Draft version February 2, 2008 

Preprint typeset using I4TgX style emulateapj v. 6/22/04 



IMPACT OF DARK MATTER SUBSTRUCTURE ON THE MATTER AND WEAK LENSING POWER SPECTRA 

Bradley Hagan 1 , Chung-Pei Ma 2 , Andrey V. Kravtsov 3 

Draft version February 2, 2008 



in 
o 
o 

(N 



in 

(N 
> 

r- : 
in . 
in - 
^t- ■ 

o : 

in . 
o ■ 

^ ' 

6 ■ 
a: 



ABSTRACT 

We explore the effect of substructure in dark matter halos on the power spectrum and bispectrum of matter 
fluctuations and weak lensing shear. By experimenting with substructure in a cosmological N = 512 3 simu- 
lation, we find that when a larger fraction of the host halo mass is in subhalos, the resulting power spectrum 
has less power at 1 < k < 100/; Mpc" 1 and more power at k > 100/i Mpc" 1 . We explain this effect using an 
analytic halo model including subhalos, which shows that the 1 < k < 100/t Mpc" 1 regime depends sensitively 
on the radial distribution of subhalo centers while the interior structure of subhalos is important at k > 100/z 
Mpc -1 . The corresponding effect due to substructures on the weak lensing power spectrum is up to ~ 1 1% at 
angular scale I < 10 4 . Predicting the nonlinear power spectrum to a few percent accuracy for future surveys 
would therefore require large cosmological simulations that also have exquisite numerical resolution to model 
accurately the survivals of dark matter subhalos in the tidal fields of their hosts. 
Subject headings: cosmology: theory — dark matter — large-scale structure of the universe 



1. INTRODUCTION 

One phenomenon to emerge from A^-body simulations of 
increasingly higher resolution is the existence of substructure 
(or subhalos) in dark matter halos (e.g., Tormen et al. 1998; 
Klypin et al. 1999b; Moore et al. 1999; Ghigna et al. 2000). 
These small subhalos, relics of hierarchical structure forma- 
tion, have accreted onto larger host halos and survived tidal 
forces. Depending on their mass, density structure, orbit, and 
accretion time, the subhalos with high central densities can 
avoid complete tidal destruction although many lose a large 
fraction of their initial mass. These small and dense dark 
matter substructures, however, are prone to numerical arti- 
facts and can be disrupted due to insufficient force and mass 
resolution. Disentangling these numerical effects from the 
actual subhalo dynamics is an essential step towards under- 
standing the composition and formation of structure. Quan- 
tifying the effects due to dark matter substructure is also im- 
portant for interpreting weak lensing surveys, which are sen- 
sitive to the clustering statistics of the overall density field. 
The level of precision for which surveys such as SNAP are 
striving (Massey et al. 2004) suggests that theoretical predic- 
tions for the weak lensing convergence power spectrum need 
to be accurate to within a few percent over a wide of range of 
scales (e.g. Huterer & Takada 2005). At this level, subhalos 
may contribute significantly to the nonlinear power spectrum 
because they typically constitute about 10% of the host mass. 

In the sections to follow, we examine the effects of sub- 
structure on the matter and weak lensing power spectra with 
two methods. In § 2 we use the result of a high resolution N- 
body simulation and quantify the changes in the power spec- 
tra when we smooth out increasing amounts of substructures. 
Our other approach, detailed in § 3, is to incorporate substruc- 
ture into the analytic halo model. The results are dependent 
on the parameters used in the model, but they provide useful 
physical insight into the results from Af -body simulations. We 
summarize and discuss the results in § 4. 
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2. SUBSTRUCTURE IN SIMULATIONS 

We use the outputs of a cosmological dark-matter-only 
simulation that contains a significant amount of substruc- 
ture. This simulation is a concordance, flat ACDM model: 
Q m = 1 - £1a = 0.3, h = 0.7 and erg = 0.9. The box size is 
120/T 1 Mpc, the number of particles is 512 3 , and the particle 
mass is 1.07 x 10 9 /i _1 Mq. The simulation uses the Adaptive 
Refinement Tree A/-body code (ART; Kravtsov et al. 1997; 
Kravtsov 1999) to achieve high force resolution in dense re- 
gions. In this particular run the volume is initially resolved 
with a 1024 3 grid, and the smallest grid cell found at the end 
of the simulation is 1.8/i _1 kpc. The actual resolution is about 
twice this value (Kravtsov et al. 1997). More details about the 
simulation can be found in Tasitsiomi et al. (2004). 

To quantify the effects of subhalos on the matter and weak 
lensing power spectra, we first identify the simulation parti- 
cles that comprise subhalos within each halo. This is achieved 
using a version of the Bound Density Maxima algorithm 
(Klypin et al. 1999a), which identifies all local density peaks 
and therefore finds both halos and subhalos. It identifies the 
particles that make up each of the peaks and removes those 
not bound to the corresponding halo. As a controlled exper- 
iment, we then smooth out the subhalos within the virial ra- 
dius of each host halo by redistributing these subhalo parti- 
cles back in the smooth component of the host halo according 
to a spherically-symmetric NFW profile (Navarro, Frenk, & 
White 1996). For the concentration parameter c of the profile, 
we do not use the fitting formulae (e.g., Bullock et al. 2001; 
Dolag et al. 2004) but instead fit each host halo individually to 
take into account the significant halo-to-halo scatter in c. We 
therefore smooth over the subhalos and increase the normal- 
ization but not the shape of the spherically averaged profile 
of the smooth component to accommodate the mass from the 
subhalo component. 

This smoothing procedure also serves as a simple model for 
the effects of resolution on the abundance of subhalos in sim- 
ulations, in which the lack of sufficient resolution will cause 
an incoming small halo to be disrupted quickly and lose most 
of its particles over its short-lived orbit. We quantify this ef- 
fect by experimenting with different cut-offs on the subhalo 
mass: subhalos with masses below the cut-off are removed 
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FIG. 1. — Effects of dark matter substructure on the matter fluctuation 
power spectrum (top) and the equilateral bispectrum (bottom) of a N = 512 3 
cosmological simulation. Plotted is the ratio of the original spectrum to that 
with a subset of the substructures smoothed out. The curves (from top down) 
correspond to increasing subhalo mass cut-offs, below which the mass in the 
subhalos is redistributed smoothly back into the host halo. The bottom panel 
is plotted to a lower k because B(k) becomes too noisy. The effects due to 
substructures are up to ~ 12% in A 2 and 24% in B. 



and have their particles spread over the host halo; subhalos 
with masses above the cut-off are left alone. Increasing mass 
cut-offs should roughly mimic increasingly lower resolution 
simulations because only higher mass subhalos will stay in- 
tact in the halo environment. 

We then calculate the matter fluctuation power spectrum 
P(k), the Fourier transform of the 2-point correlation function, 
and the matter bispectrum B{ki,k2,kj,), the Fourier transform 
of the 3-point correlation function. In order to compute P 
and B at large k without using an enormous amount of mem- 
ory, we subdivide the simulation cube into smaller cubes and 
stack these on top of each other (called "chaining the power" 
in Smith et al. 2003). A typical stacking level used is 8, mean- 
ing that we subdivide the box into 8 3 cubes and stack these. 
We use stacked spectra for the high-A regime and unstacked 
spectra for low-k. Finally, we subtract shot noise (oc l/N) 
from the outputted spectra to eliminate discreteness effects. 

Fig- n shows the effects of substructure on the dimension- 
less power spectrum A 2 (k) = A 3 P(k) / (2ir 2 ) and the equilat- 
eral bispectrum B{k\) {k\ = ki = A13). Plotted is the ratio of 
the spectrum from the raw ART output divided by the spec- 
trum from the altered data. The altered data have no subha- 
los with masses below the labeled mass cutoff. The mass in 
the removed subhalos has been redistributed smoothly into the 



host halo as described above. The deviation from the original 
spectrum becomes larger as the cutoff is increased because 
more subhalos have been smoothed out. For a given cutoff, 
the figure shows that a simulation with dark matter substruc- 
tures (such as the raw ART output) has more power at k > 100 
h Mpc" 1 and less power at 1 < k < 100 h Mpc" 1 than a simu- 
lation with smoother halos. We believe these opposite behav- 
iors reflect the two competing factors present in our numeri- 
cal experiments: removal of mass within subhalos, which af- 
fects scales comparable to or below subhalo radii (and hence 
k > 100/z Mpc" 1 ), and addition of this mass back into the 
smooth component of the halo, which affects the larger scales 
of 1 < k < 1 00 h Mpc" 1 . The ratio approaches unity for k<lh 
Mpc" 1 simply because the mass distribution on scales above 
individual host halos is unaltered. We will examine these ef- 
fects further in the context of the halo model in § 3. 

For a given curve in Fig. ^ we nave a l so calculated the 
contributions from subhalos in host halos of varying masses to 
quantify the relative importance of cluster versus galactic host 
halos. For the 10 12 ■ 5 /i" 1 M curve, e.g., we find that smooth- 
ing over the subhalos in host halos above 10 14 and 10 13 M Q 
account for 5% and 10% in the total 12% dip seen in Fig.^ 
respectively. For the 10 11 s h~ x MQ cutoff, the numbers are 2% 
and 5% of the total 6% dip. 

The halos found in A^-body simulations are generally triax- 
ial. When we redistribute the subhalo particles, however, we 
assume for simplicity a spherical distribution. This assump- 
tion makes the altered halos slightly rounder. One can esti- 
mate how this effect changes the power spectrum by using the 
halo model without substructure. Smith & Watts (2005) in- 
corporated a distribution of halo shapes found by Jing & Suto 
(2002) from cosmological simulations into the halo model (ig- 
noring the substructure contribution). Compared with the case 
where all the halos are spherical, they observed a peak decre- 
ment in the power spectrum of about 4% for Awl Mpc" 1 . 
The corresponding effect in our calculations would be much 
smaller since we redistribute only the subset of particles that 
belong to subhalos into the rounder shape (e.g., about 10% of 
all particles in the case where the cutoff was 10 12 5 /! _1 Mq). 
Thus, by extending the results of Smith & Watts, we expect 
the spurious rounder halos in our study to account for less 
than 0.5% of the total 12% drop. 

Fig-Elshows the weak lensing convergence power spectrum 
A 2 (Z) corresponding to the matter power spectrum A 2 (k) in 
Fig-E It is calculated from A 2 (k) using Limber's approxima- 
tion and assumption of a flat universe: 

kW M H c 



X 3 d X ^MA 2 (k = e/ X ,a). 



a (X) 

(1) 

Here X lS the comoving radial distance, and the weak lens- 
ing weight is W( X ) = f^"" dx'p(x')(x' ~X)/x', where P(X) 
is the distribution of source galaxies such that J p(x)dx = 1 
(see, e.g., Bartelmann & Schneider 2001). Here we assume 
for simplicity that all sources are at one redshift (z = 1) and 
use the z = simulation output to estimate A 2 (k). Host ha- 
los at higher redshift may have a larger fraction of their mass 
in substructures because the typical subhalo accretion epoch 
would be more recent and there would be less time for tidal 
disruption. The effects of z ~ 0.5 substructure may therefore 
be somewhat larger than shown here for z = 0, although we do 
not expect the change to be significant. 

The gray band in Fig. [2] marks the extent of uncertainties 
from sample variance and shape noise in weak lensing mea- 
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FIG. 2. — Effects of substructure on the weak lensing convergence power 
spectrum for the same subhalo mass cutoffs as in Fig. 1. The sources are 
assumed to be at redshift 1 . The gray band is the 1-cr statistical error assuming 
Gaussian fields. We take f s ty — 0.25, 7™j = 0.2, n= 100/arcmkr, and a band 
of width £/l0. 



surements assuming Gaussian density fields (Kaiser 1998): 
a(Al) I 2 



(2£+ l)f sky V 2™A 2 K 



i + 



(2) 



where f s ^ is the fraction of the sky surveyed, 7, mj is the rms 
ellipticity of galaxies, and h is the number density of galaxies 
on the sky. The error is dominated by the sample variance 
on large scales (first term in eq. U) and the "shape noise" on 
small scales. Our assumption of Gaussianity is not applicable 
for the angular scales shown in the plot because the scales 
plotted are near or below the size of individual halos, but the 
errors shown should be a useful reference and have been used 
in previous studies. A reliable estimate of the error would 
presumably require a ray tracing calculation which is beyond 
the scope of this paper. 

3. SUBSTRUCTURE IN THE HALO MODEL 

To gain a deeper understanding of the simulation results in 
Fig s -E an d|2] we use the semi-analytic halo model to build 
up the nonlinear power spectrum from different kinds of pairs 
of mass elements that may occur in halos (e.g. Ma & Fry 
2000; Peacock & Smith 2000; Seljak 2000; Scoccimarro et 
al 2001). The original halo model assumes that all mass re- 
sides in virialized, spherical halos without substructures. One 
can then build the matter power spectrum from the different 
kinds of pairs of particles that contribute to the 2-point clus- 
tering statistics by writing P(k) = P\h{k) + P2h{k), where the 1- 
halo term P\h contains contributions from particle pairs where 
both particles reside in the same halo, and the 2-halo term Py, 
is from pairs where the two particles reside in different halos. 
The 1-halo term is a mass-weighted average of single halo 
profiles and dominates on the scales of interest (k > lh Mpc" 1 ) 
in Fig.^because close pairs of particles are more likely to be 
found in the same halo. The 2-halo term is closely related to 
the linear power spectrum and is important only at large sepa- 
ration (i.e. small k) where a pair of particles is more likely to 



be found in two distinct halos. Similarly, the bispectrum can 
be constructed from the different classes of triplets of particles 
(see, e.g., Ma & Fry 2000). 

The original halo model can be readily extended to take into 
account a clumpy subhalo component in an otherwise smooth 
host halo. Sheth & Jain (2003), e.g., decompose the original 
1-halo term into Pu, = P ss + P sc + P\ c + Pic, where "s" denotes 
smooth and "c" denotes clump. The smooth-smooth term, P ss , 
arises from pairs of particles that both belong to the smooth 
component of the same host halo. This term is identical to the 
original 1-halo term except for an overall decrease in ampli- 
tude by the factor (1 -/) 2 , where / is the fraction of the total 
halo mass that resides in subhalos. The smooth-clump term, 
P sc , is due to having one particle in a subhalo (clump) and the 
other in the host halo (smooth). The 1- and 2-clump terms, P\ c 
and Pic, come from having both particles in the same subhalo 
and in two different subhalos, respectively. Explicitly, 

I ' dM N(M)M 2 U 2 (k,M) (3) 



P ss (k) : 



Psc(k) = 2 



Pidk) ■ 



P 

(1-/) 



dM N(M)MU(k,M)U c (k,M) 



dm n(m,M)mu(k,m) 



(4) 



1 



dMN(M) / dm n(m,M)m u (k,m) (5) 



(6) 



Pic(k)= ~z I dMN(M)U c z (k,M) 



dmn(m,M)mu(k, m) 



where U(k,M), u(k,m), and U c (k,M) are the Fourier trans- 
forms of the host halo radial density profile, the subhalo radial 
density profile, and the radial distribution of subhalo centers, 
respectively. N(M)dM gives the number density of host ha- 
los with mass M, and n{m,M)dm gives the number density of 
subhalos of mass m inside a host halo of mass M. A similar 
expression can be written down for the 2-halo term P21,, which 
is also included in our calculations. 

Fig-E]illustrates the contributions from the individual terms 
in the halo model. We use the NFW profile (truncated at 
the virial radius) for the input host halo U(k,M) and subhalo 
u(k,m), and the concentration c(M) = co(M/lO 14 M )"° 1 with 
co = 1 1 that we find to approximate the ART host halos and 
is identical to Dolag et al (2004) except for a 15% increase 
in amplitude. We use Cg ub = 3 for the subhalos to take into 
account tidal stripping but also compare different values in 
Fig. 4 below. For the distribution of subhalo centers, U c (k), 
we compare the profile of NFW with that of Gao et al. (2004), 
who find the number of subhalos within a host halo's virial ra- 
dius r v to be 

N«x) (l + ac) X P 

= , x=r/r v (7) 

Ntot l+acx a ' ' 

where a = 0.244, a = 2,(3 = 2.75, c = r v /r s , and N, or is the to- 
tal number of subhalos in the host. Since this distribution at 
small r is shallower (oc r~~° ,25 ) than the inner part of the NFW 
profile (oc r" 1 ), its Fourier transform U c (k) at high k is about 
a factor of 10 lower than that of the NFW profile. This decre- 
ment results in a much lower A sc and Ai c as shown in Fig. [5] 
(dashed vs dotted curves). We find the subhalo centers in the 
ART simulation to follow approximately the distribution of 
Gao et al. although there is a large scatter. 
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FIG. 3. — Comparison of individual subhalo terms in the halo model. 
The smooth-clump A 2 C and 2-clump A3, terms depend on the distribution 
of subhalo centers U c within a host halo, having a much lower amplitude at 
k > Ih Mpc -1 when U c has the cored isothermal profile of Gao et al. (2004) 
compared with the cuspy NFW profile. The smooth-smooth and 1-clump 
A| terms are independent of U c - See text for parameters used in the model. 



We use the mass function of Sheth & Tormen (1999) for 
the host halos N(M) and a power law n(m,M) oc m~ that 
well approximates the subhalo mass function in the ART sim- 
ulation. The latter is normalized so that the total mass of 
subhalos in a host halo adds up to / times the host mass M 
(/ = 0.14 in Fig. 3). To compare the halo model with simula- 
tions, we set the lower limit on the P ss integral to 10 10 /i _1 Mq, 
which is the smallest halo present in the simulation (about ten 
times the simulation particle mass). The lower limit on the 
outer integrals of P sc , P\ c , and P^ c corresponds to the smallest 
halo that contains substructure, which we set to the small- 
est halo that we considered for erasing substructure in § 2: 
2 x 10 12 /r' Mq. Similarly, the lower limit on the inner inte- 
grals of these terms is set to the smallest subhalo that can be 
resolved (i.e. I0 m h~ l Mq). 

Fig. 4 compares the sum of all the terms in the halo model 
with the simulation result from Fig. 1. As in Fig.^ we il- 
lustrate the effects due to substructures by dividing out the 
power spectrum from the original (smooth) halo model, i.e., 
Smooth = £ 3 A s /[(1-/) 2 (2tt 2 )], where P ss is given before in 
eq. (|4). Fig.[4]shows that the halo model is able to reproduce 
qualitatively the simulation results when the subhalo centers 
in the halo model are assigned the shallower distribution of 
Gao et al. The feature of A 2 sub / A 2 mwoth < 1 at 1 < * < 100 h 
Mpc" 1 is mainly caused by the drop in the smooth-clump term 
relative to the smooth-smooth term at k > 1 h Mpc" 1 shown in 
Fig.|3] The ratio A 2 ulh / A 2 smooth becomes > 1 only at k > 100 h 
Mpc" 1 when the 1-clump term finally takes over. Fitting the 
halo model to actual simulation results is clearly not exact in 
part due to the large scatters in the properties of simulated ha- 
los, e.g., the concentration (for both hosts and subhalos), the 
subhalo mass fraction /, and the maximum subhalo mass in 
each host halo. The halo model allows us to study the depen- 
dence of clustering statistics on these parameters (see Fig. 4). 
In addition, a number of effects are neglected in the current 
halo model, e.g., tidal effects are likely to reduce the number 
of subhalos (modeled by U c of Gao et al. here) as well as their 
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FIG. 4. — Same ratio of A as in Fig. Q] but comparing simulation 
(symbols; same 10 12 5 /j~'Mq curve in Fig. QJ with halo model predictions 
(plain curves). The two agree qualitatively when the shallower distribu- 
tion of Gao et al. for subhalo centers U c (k) is used (bottom 3 curves) in 
the halo model (but not for an NFW U c (k); dotted). The detailed model 
prediction depends on halo parameters: the solid curve uses the same pa- 
rameters as in Fig. 3; the dashed shows how a larger subhalo concentration 
(c sub = qJ ,b (M/lO 14 M r° 1 ; c™ b = 1 1 vs 3) steepens the curve at high k; the 
dash-dotted shows how a smaller subhalo mass fraction / (0.1 vs. 0.14) raises 
the dip. 



outer radii (not modeled here) towards host halo centers; a 
larger amount of stripped subhalo mass may also be deposited 
to the inner parts of the hosts, resulting in a radius-dependent 
subhalo mass fraction / within the host. 

Fig. 4 also shows that the halo model predicts the oppo- 
site effect due to substructure (i.e. A 2 uh j A 2 mooth > 1 at all k) 
if the subhalo centers U c (k) are assumed to follow the NFW 
distribution like the underlying dark matter. The sign of this 
effect is consistent with the previous subhalo model study of 
Dolney et al. (2004), which assumed the same NFW profile 
for the subhalo centers and the hosts and obtained a matter 
power spectrum that had a higher amplitude for all k when the 
substructure terms were included. Their results differ slightly 
from ours because of different integration limits. Subhalos 
in recent simulations like that of Gao et al., however, show 
a much shallower radial distribution in the central regions 
of the host halos, and inclusion of gas dynamics appears to 
have little effect on the survivability of subhalos (Nagai & 
Kravtsov 2005). The shallow distribution is apparently due to 
tidal disruptions, even though the precise shape of the distri- 
bution is still a matter of debate (e.g., Zentner et al. 2005). 
We have also experimented with a third distribution U c (k) that 
has the NFW form but is less concentrated. We are able to 
bring A 2 ub / A 2 mooth below unity only when the concentration 
is reduced by a factor of more than 2.5, and only when this 
reduction factor is increased to ~ 100 would we get a com- 
parable dip as the curves for ART simulation and Gao et al. 
in Fig.|4] It is interesting to see if we can mimic the behav- 
ior of the ART simulation without using subhalos in the halo 
model. We try replacing the one-halo term, Pi/,, by one that 
is a simple superposition of a Gao et al. profile and the usual 
NFW profile. This accounts for the fact that ~ 90% of the 
mass is in a smooth NFW profile and that ~ 10% is in sub- 
halos, which follow a flatter profile. One would not expect 
the high-fe regime to agree as it is dominated by the subha- 
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los (the 1 -clump term, specifically). The intermediate range, 
1 < k < 100/; Mpc -1 , is dominated by the host halo itself, but 
we find no similarity in this range either. Subhalos are there- 
fore needed if the halo model is to recreate the ART results. 

4. DISCUSSION 

The purpose of this work is to provide a physical under- 
standing of the effects of substructures on clustering statistics. 
By experimenting with dark matter substructures in a cosmo- 
logical simulation with 51 2 3 particles, we have shown that the 
power spectra of matter fluctuations and weak lensing shear 
can change by up to ~ 12% (and up to ~ 24% in the bispec- 
trum) if a significant amount of substructures is not resolved 
in a simulation. When a larger mass fraction of the host halos 
is in the form of lumpy subhalos, we find the effect is to lower 
the amplitude of the matter and weak lensing power spectra 
at the observationally relevant ranges of k ~ 1 to 100 h Mpc -1 
and I < 10 5 , and to raise the amplitude on smaller scales. A 
similar drop in power is also seen in our analytic halo-subhalo 
model when the subhalo centers U c within a host halo are dis- 
tributed with a shallower radial profile than the underlying 
dark matter (as expected due to tidal effects). A way to un- 
derstand the drop involves looking at where the dense regions 
are. When U c has an NFW form the subhalos basically trace 
the smooth background. Thus, there is never a decrease in 
power when the smooth-smooth and smooth-clump terms are 
added because dense regions are in nearly the same relative 
positions. When we use a shallower profile for U c , the sub- 
halos are not as numerous in the denser inner regions of the 
background halo. This decrease in the overlap between dense 
clump regions and the dense inner regions causes the drop in 
power. 

We have quantified the effects of substructures on clustering 
statistics by erasing substructures in an N = 512 3 simulation. 
An important related question is whether N = 512 3 , single- 
mass resolution simulations such as the one used in our study 
has sufficient resolution to measure the power spectrum to the 
few-percent accuracy required by future surveys. Note that 
at least hundreds of particles and force resolution of ~ kilo- 
parsec are required to ensure subhalo survival against tidal 
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forces, placing stringent requirements on the dynamic range 
of simulations. Multi-mass resolution simulations designed 
for subhalo studies, on the other hand, do not give reliable 
predictions for P(k) on quasi-linear scales due to compro- 
mised resolution outside highly clustered regions. The fact 
that the curves in Figs. 1 and 2 continue to change at the 
few-percent level each time the mass threshold is lowered by 
0.5 dex from 10 12 ' 5 to 1O 1O - 5 A _1 M suggests that subhalos of 
M < 1O 1O ' 5 /i _1 M0 may still be affecting the power spectra at 
a comparable level and that N > 512 3 would be required. We 
also find > 3% changes in P(k) at k ~ \0h Mpc" 1 in the halo 
model as the minimum subhalo mass in the integration limit 
of eq. (3) is lowered to 10 7 h~ l M Q (although the exact predic- 
tions are sensitive to the slope of the subhalo mass function, 
which is assumed to be -1 .9 here.) Careful convergence stud- 
ies with higher resolution aided by insight from this study and 
detailed semi-analytic models for halo substructure will likely 
be needed to determine N. 

There are other challenges to predicting accurately the weak 
lensing signal on single halo scales. The effect of neutrino 
clustering could cause a rise in weak lensing convergence of 
~ 1% at £ ~ 2000 (Abazajian et al. 2005). Two recent groups 
have investigated different aspects of baryon effects. White 
(2004) found that baryonic contraction and its subsequent im- 
pact on the dark matter distribution is capable of causing an 
increase in the weak lensing convergence power of a few per- 
cent at i > 3000. Zhan & Knox (2004), on the other hand, use 
the fact that the hot intracluster medium does not follow the 
dark matter precisely and predict an opposite effect: a sup- 
pression of weak lensing power of a few percent at £ > 1000. 
Unlike the effects of substructure and neutrino clustering, the 
baryon effects cause departures from the pure dark matter 
weak lensing signal that only get larger with increasing I. 
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